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ABSTRACT 


This  report  is  concerned  with  the  design  of  a  bounded 
wave  Electromagnetic  Pulse  (EMP)  simulator.  Different  types 
of  simulators  are  described  and  their  pros  and  cons  are 
discussed.  A  detailed  design  is  given  for  a  wire  grid 
parallel  plate  type  simulator.  The  fields  inside  the 
simulator  as  well  as  the  radiated  fields  around  the  simulator 
are  also  computed. 


RESUME 


Ce  rapport  traite  de  la  conception  d'un  simulateur 
d' impulsions  electro-magnetiques  ( IEM )  de  type  h  onde 
guidee.  Differents  types  de  simulateurs  sont  revus,  ainsi 
que  leurs  avantages  et  leurs  inconvenients .  Les  plans 
detailles  d'un  simulateur  de  type  i  grille  parallels  y  sont 
decris.  Les  champs  electro-magnetiques  a  l'interieur  et  a 
l'exterieur  du  simulateur  sont  aussi  calcules. 
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EXECUTIVE  SUMMARY 


Electromagnetic  Pulse  (EMP)  simulators  are  used  to 
simulate  the  EMP  generated  by  a  nuclear  weapon  and  to  harden 
equipment  against  the  effects  of  EMP. 


At  present,  DREO  has  aim  EMP  simulator  for  testing 
computer  terminal  size  equipment.  In  order  to  develop  the  R&D 
capability  for  testing  larger  objects  such  as  a  helicopter, 
a  much  bigger  threat  level  facility  is  required. 


This  report  concerns  the  design  of  a  bounded  wave  EMP 
simulator  suitable  for  testing  large  size  equipment. 
Different  types  of  simulators  are  described  and  their  pros  and 
cons  are  discussed.  A  bounded  wave  parallel  plate  type 
simulator  is  chosen  for  it's  efficiency  and  the  least 
environmental  impact. 


Detailed  designs  are  given  for  6  m  and  10  m  parallel 
plate  type  wire  grid  simulators.  Electromagnetics  fields 
inside  and  outside  the  simulators  are  computed.  Preliminary 
specifications  for  a  pulse  generator  required  for  the 
simulator  are  also  given.  Finally,  the  electromagnetic  fields 
radiated  from  the  simulator  are  computed  and  discussed. 
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DESIGN  OF  A  BOUNDED  WAVE  EMP  SIMULATOR 


Intended  as  second  stage  simulator  for  DREO 


1.0  INTRODUCTION 

The  Phase  I  simulator  in  the  EMP  Project  of  Defence  Research 
Establishment  Ottawa,  started  in  October  1987,  is  installed  and  is  now  being 
tested.  It  is  a  small  laboratory  simulator  (height  1  metre),  installed  in 
building  39  and  intended  for  experiments  which  are  expected  to  continue 
through  FY  89/90. 

In  Phase  II  a  larger  EMP  simulator  with  a  larger  test  volume  is  needed. 
This  second  stage  simulator,  intended  as  a  R&D  tool  for  EMP  Protection 
Studies,  will  be  built  in  the  field  behind  building  39.  Some  of  the 
components  of  the  Phase  I  simulator  will  be  used  for  the  second  stage 
simulator,  such  as  the  shielded  room,  the  sensors  and  the  data  acquisition  and 
processing  instruments.  New  equipment  will  include  the  criterion  level  pulse 
generator,  the  transmission  line  with  its  support  structure  and  a  fibre-optic 
data  transmission  system.  If  possible  the  100  kV  pulser  from  the  small 
simulator  will  also  be  used  for  low  level  survey  studies  until  the  larger 
pulser  is  available. 

Previous  DREO  proposals  for  a  second  stage  simulator  specified  a  6  m 
bounded-wave  simulator.  This  height  is  sufficient  to  accept  vehicle-sized 
systems.  Larger  bounded-wave  simulators  were  considered  too  expensive. 

Pree-field  simulators  are  expected  to  have  unacceptable  radiation  for 
the  proposed  site.  Therefore  a  choice  has  to  be  made  among  the  different 
types  of  bounded-wave  simulators.  The  dimensions  should  be  limited  to  an 
acceptable  size  and  the  shape  adapted  to  the  available  site. 

In  this  report  the  general  rules  for  designing  a  bounded-wave  simulator 
are  given,  with  two  examples,  one  with  a  height  of  6m  and  another  with  a 
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height  of  10m,  respectively  called  a  6m  and  10m  line.  The  calculations  of  the 
expected  radiation  are  given  in  Appendix  A. 

The  work  has  been  performed  as  a  part  of  the  Consulting  Service 
Contract  mentioned  in  the  agreement  attached  to  letter  DREO  3743A-GE  (ED)  of 
9  May  1988. 
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2.0  SIMULATOR  REQUIREMENTS 


As  starting  point  for  the  design  the  following  restrictions  are  made: 


-  Only  the  exo-atmospheric  EMP  has  to  be  simulated. 

-  The  time  history  of  this  pulse  should  be  as  described  in  appropriate 
NATO  documents. 

-  The  working  volume  should  have  a  length  of  about  20m,  a  maximum  width 
of  20m  and  a  height  not  more  than  10m. 

-  If  possible  the  dimensions  of  the  working  volume  should  be  flexible. 

-  The  maximum  field  strength  should  be  50-65  kV/m. 

-  The  radiation  should  be  as  low  as  possible. 
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3.0  SELECTION  OF  SIMULATOR  TYPE 

Since  the  available  site  does  not  permit  a  high  radiation  level  the 
bounded -wave  (or  transmission  line)  EMP  simulator  is  the  most  suitable  type 
for  this  purpose.  An  advantage  is  that  this  type  of  simulator  can  also  be 
used  for  endo-atmospheric  EMP  simulation,  perhaps  important  in  a  later  stage. 

It  should  be  kept  in  mind  that  it  is  impossible  to  completely  duplicate 
the  EMP  threat  fields  by  simulation.  Each  type  of  simulator  has  its 
limitations.  For  the  bounded-wave  simulator  they  are: 

-  vertical  polarization  and  horizontal  propagation 

-  no  ground  interactions  (a  limitation  for  ground-based  systems,  an 
advantage  for  air  born  systems) 

-  potential  interaction  of  simulator  and  object  under  test 

-  limited  working  volume 

Special  features  are  the  high  threat  level  at  relatively  low  cost  and 
the  proper  wave  impedance.  EMP  simulation  with  horizontal  polarization  has 
been  left  out  of  consideration  here,  because  of  the  high  cost. 

A  selection  can  be  made  among  the  following  types  of  bounded-wave 
simulators,  see  fig.  1: 

-  parallel  plate  simulator  (PPS) 

-  conical  plate  simulator  (CPS) 

-  triangular  plate  simulator  (TPS) 

-  rhombic  wire  simulator  (RVS) 


3.1  Parallel  Plate  Simulator 


This  type  of  simulator  operates  by  guiding  an  electromagnetic  vave 
across  «  test  object  situated  between  tvo  metallic  surfaces  of  a  transmission 
line,  operating  in  the  transverse  electromagnetic  (TEH)  propagation  mode.  That 
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is  a  traveling  wave  with  a  vertical  electric  field  and  a  horizontal  magnetic 
field.  The  transmission  line  is  driven  at  one  end  by  a  generator  and  loaded 
at  the  other  end  by  a  termination.  The  generator  and  termination  are 
connected  to  the  horizontal  transmission  line  section  by  transition  sections 
(tapers).  The  characteristic  impedance  of  the  transmission  line  is  determined 
by  the  width-to-height  ratio,  a/h,  and  can  be  made  constant  from  generator  to 
termination  by  tapering  the  transitions.  When  the  transmission  line  is 
correctly  terminated,  R  =  ZQ,  all  energy  delivered  by  the  generator  should  be 
dissipated  by  the  terminating  resistor  R.  Without  reflections  there  will  be 
only  traveling  waves  and  no  standing  waves. 

However  the  TEM  mode,  launched  in  the  input  taper,  experiences 
reflections  before  reaching  the  termination.  At  the  junction  of  the  tapers 
with  the  horizontal  section,  energy  is  converted  to  other  propagating  modes. 
Moreover,  for  frequencies  higher  than  the  frequency  at  which  plate  spacing  h 
is  a  half  wavelength,  it  is  possible  to  maintain  higher  order  modes  in  the 
transmission  line.  In  the  literature  deep  notches  in  the  frequency  spectrum 
of  the  field  in  the  working  volume  are  reported  (2  to  9). 


3.2  Conical  Plate  Simulator 


To  prevent  reflections  caused  by  local  deviations  of  ZQ  at  the 
junctions  of  the  horizontal  section  of  the  PPS  and  to  avoid  radiation 
generated  in  this  section  the  idea  was  born  to  eliminate  this  horizontal  part 
of  the  simulator  and  to  use  the  input  taper  as  working  volume.  The  question 
is  then  how  to  prevent  reflections  at  the  end  of  the  input  taper.  The  second 
step  was  to  eliminate  the  output  taper  as  well  and  to  leave  the  end  of  the 
generator  taper  open.  The  wave  launched  by  the  generator  in  the  taper  should 
then  be  radiated  in  free  space.  In  this  case  the  taper  behaves  like  a  horn 
antenna.  Since  the  dimensions  of  the  horn  antenna  are  too  small  for  radiating 
low  frequencies,  some  kind  of  termination  is  needed  at  the  end.  This  type  of 
simulator  is  called  conical  plate  simulator  (CPS). 
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3 • 3  Triangular  Plate  Simulator 

In  rhe  TPS  only  the  horizontal  parallel  plate  section  of  the  PPS  is 
eliminated.  The  input  taper  is  directly  connected  to  the  output  taper  in  the 
middle  of  the  simulator.  It  is  obvious  that  the  discontinuity  between  the 
tapers  is  larger  than  in  the  PPS,  with  larger  reflections  as  result. 


3 . 4  Rhombic  Wire  Simulator 

Since  a  major  part  of  the  current  is  concentrated  near  the  edges  of  the 
plates  King  and  others  replaced  the  tapers  of  the  TPS  by  a  wire  or  tube  at  the 
two  edges  of  the  tapers,  thus  forming  a  rhombic  antenna  above  the  ground  (4  to 
7).  This  structure  is  mainly  used  for  R&D  purposes  since  the  radiated 
electromagnetic  fields  of  such  antennae  can  accurately  be  determined.  Uith 
this  model  the  analysis  of  phenomena  can  be  used  for  other  EMP  simulators. 


3.5  Examination  of  the  Pros  and  Cons 


Let  us  go  back  to  the  PPS.  The  reflections,  standing  wave  ratios 
(SVR's)  and  notches  reported  in  the  literature  were  calculated  and  measured 
with  continuous  wave  (CV)  excitation  (2,3).  CV  excitation,  however,  is 
different  from  a  single  shot,  pulse  excitation.  Vith  CU  excitation  the 
discontinuities  experienced  by  the  traveling  waves  cause  reflections  and  since 
the  excitation  is  continuous,  standing  waves  are  formed.  The  SVR  is  dependent 
on  the  mismatch. 

Vith  a  single  pulse  excitation  no  standing  waves  will  exist  and  there 
is  no  SVR.  Radiation  of  the  resonant  structure,  existing  with  CV  exitation  at 
frequencies  related  to  the  simulator  height  h,  causing  notches  inside  the 
wavegjlde,  cannot  build  up  large  magnitudes  with  a  single  pulse  excitation. 

For  pulse  excitation  the  time-space  concept  can  better  be  used  for 
describing  the  electromagnetic  properties  of  EMP  Simulators  (See  fig.  2).  The 
clear-time,  defined  as  the  time  from  the  initial  observation  of  the  wave  to 
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the  time  of  arrival  of  the  first  interference  wave  (reflection)  at  the  point 
of  observation,  is  a  more  useful  parameter.  It  may  be  that  at  frequencies 
where  notches  are  observed  with  CV  excitation,  some  energy  will  be  radiated 
with  pulse  excitation,  however  the  guided  wave  travelling  between  the  parallel 
plates  will  not  be  influenced  by  that.  The  reflections  at  the  discontinuities 
of  the  PPS  can  be  reduced  by  designing  relatively  long  tapers  (low  elevating 
angle)  and  by  "smoothing"  the  angle  of  the  junctions. 

A  small  elevation  angle  has  more  advantages.  One  of  them  is  that  the 
wave  front  launched  in  the  taper  is  less  spherical,  causing  less  distortion  of 
the  field  orientation  near  the  top  plate  of  the  horizontal  section.  The 
spherical  wave  in  the  transistion  section  tries  to  transform  into  a  plane  wave 
in  the  working  volume  to  establish  the  TEM  propagating  wave  in  the  horizontal 
wave  guide.  For  that  reason  it  is  recommendable  to  design  a  minimum  length  of 
this  section  equal  to  3h.  With  a  good  design  the  reflections  at  the 
terminating  resistor  of  the  PPS  can  be  made  very  small  as  has  been  shown  in 
the  past  (10). 

The  advantages  claimed  for  the  CPS  are: 

-  no  reflections 

-  no  typical  radiation  frequencies 

-  no  distortion  of  the  field  orientation 

Since  the  horn  antenna  is  not  ideally  matched  to  the  impedance  of  free 
space  (ZQ  =  120n  S)  -  the  low  frequencies  cannot  be  radiated  and  have  to  be 

terminated  -  and  even  when  this  large  terminating  network  is  made  very  sparce 

will  be  seen  by  the  launched  wavefront  -  reflections  will  exist.  These 
reflections  arrive  in  the  working  volume  with  a  delay  related  to  the  clear 

time.  For  frequencies  with  wavelength  comparable  to  the  height  of  the  large 

terminating  network,  resonances  exist  with  consequent  resonant  radiation.  This 
can  be  reduced  by  decreasing  the  inductance  of  the  network  (11).  Disadvantages 
of  the  CPS  compared  to  the  PPS  are  the  larger  radiation  and  the  geometrical 
1/R  attenuation  of  the  launched  wave  amplitude  with  distance  from  the 
generator. 
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Comparing  the  characteristics  of  the  PPS  and  CPS  it  can  be  concluded 
that  the  PPS  is  the  best  choice  as  second  stage  EMP  simulator  for  DREO. 


I 
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4.0  THREAT  ENVIRONMENT  (17) 

A  commonly  accepted  worst-case  model  for  the  exo-atmospheric  EMP  is: 

E(t)  =  Eq  (e~at  -  e"0t)  (1) 

vi  th: 

E  =  5.3  •  104  V/m 
o 

a  =  3.7  •  106  s"1 
0  =  3.9  •  108  s"1 

The  spectrum  of  (1)  can  be  calculated  with: 

E(w)  =  E  f  -4-  -  ^4-  )  (2) 

'  o  V  a+3<o  0*-j«  ) 

With  u  =*  0  the  constant  spectral  amplitude  at  low  frequencies  and  with 
w  =  a  and  u  =  0,  the  -3  dB  cut-off  frequencies  can  be  calculated.  Figs.  3  and 
4  show  the  time  history  and  amplitude  spectrum.  The  rise  time  t  (10%  -  90%) 
=  5  ns,  the  half  width  time  t^  (50X  -  50%)  ®  200  ns,  the  peak  field  strength 
is  50  kV/m.  The  incident  field  at  the  earth  surface  can  be  considered  as  a 
plane  wave: 

E/H  -  120n  O  (3) 


The  polarization  north  and  south  of  ground  zero  (GZ)  is  horizontal  and 
at  a  large  distance  east  or  west  from  GZ  almost  vertical. 

The  analytical  model  presented  in  (1)  has  a  discontinuity  at  t  »  0  that 
physically  cannot  exist  and  that  the  EMP  will  not  have.  Moreover  this 
phenomenon  cannot  be  simulated  too. 

Therefore  the  following  expression  can  better  be  used: 
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with: 

A  =  5  •  104  V/m  .  E  .  at  t  =  0 

pk 

a  =  7.8  •  108  s"1 

8  =  3.45  •  106  s*1 

The  rise  time,  half  width  time  and  Ep^  are  the  same  as  in  model  (1), 
however  the  field  has  a  more  gradual  rise.  The  dashed  lines  in  fig.  3  and  4 
show  the  difference.  This  time  history  looks  more  like  that  of  the  EHP 
simulator.  The  relaxed  high  frequency  end  of  the  spectrum  has  important 
consequences  for  the  exposed  equipment  in  the  simulator  because  most  of  the 
responses  are  related  to  the  first  derivative  of  the  field.  This  effect  can 
be  shown  in  fig.  5,  the  first  derivative  of  the  spectrum.  Fig.  6  shows  the 
cumulative  energy  density  of  (1). 
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5.0  ELECTROMAGNETIC  FIELD  IN  THE  SIMULATOR 

The  field  in  the  parallel-plate  waveguide  of  the  simulator  is  composed 
of  TEM,  TE  and  TM  modes.  If  the  elevation  angle  0  of  the  input  taper  is 
small,  the  dominant  mode  is  a  TEM  mode.  This  mode  of  propagation  is  indistin¬ 
guishable  from  free  space  propagation. 


5.1  Field  Magnitude 

The  field  strength  in  the  simulator  is  given  by  the  expressions: 

E2  -  l  V/.  (5) 

"x  ‘  ST  *'«  <‘> 

o 

where  V  is  the  peak  pulse  amplitude,  b  the  plate  separation  in  metres 
and  Zq  the  wave  impedance  of  free  space. 

5.2  Field  Distortion 


The  dimensions  of  the  transmission  line  are  determined  by  the  size  of 
the  objects  placed  in  the  vorking  volume.  The  vertical  electric  field 
component  is  short-circuited  by  the  vertical  conductive  parts  of  the  object. 
The  height  of  the  parallel  plate  section  should  be  chosen  such  that  the 
largest  object  is  not  unduly  overstressed  by  the  increased  field  strength.  For 
field  changes  of  202  or  less  a  simulator  height  of  5  to  10  times  the  largest 
object  is  suggested.  In  practical  cases,  with  large  objects,  heights  between 
2  and  3  times  are  satisfactory.  The  transverse  dimensions  will  be  related  to 
the  number  and  width  of  the  objects  and  to  the  tolerable  field  variations 
across  the  vorking  area.  Conductors  exposed  to  the  simulated  field  reradiate 
their  induced  energy.  The  major  field  component  is  the  quasi  electrostatic 
field,  that  falls  off  inversely  with  D  .  D  is  the  distance  from  the  radiating 
conductor.  At  distances  equal  to  the  height  of  the  object,  the  field 
distortion  is  ca  502  and  at  2,  3  and  4  times  the  object  height  respectively 
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6,  2  and  IX. 


5.3  Field  Orientation 

The  electromagnetic  field  orientation  and  relative  magnitude  for  a  TEM 
wave  in  a  parallel-plate  wave  guide  can  be  displayed  by  a  plot  of  equi- 
potential  lines  and  field  lines  as  shown  in  fig.  7.  The  equipotential  lines, 
given  by  constant  v,  show  the  magnetic  field  orientation,  their  spacing 
relates  to  the  electric  field  magnitude  (volts  per  metre).  The  field  lines, 
given  by  constant  u,  show  the  electric  field  orientation,  their  spacing 
relates  to  the  magnetic  field  magnitude  (amperes  per  metre). 


5.4  Field  Uniformity 

The  field  uniformity  for  the  TEM  wave  in  the  transmission  line  depends 
or.  the  ratio  of  b/a  and  the  location  within  the  line  (2b  is  the  separation  and 
2a  the  width  of  a  parallel  two  plate  transmission  line,  see  fig.  8a).  In  fig. 
9  through  20  the  field  and  potential  distribution  is  given  for  b/a  ranging 
from  0.2  to  5.  Since  the  plots  are  symmetric,  reflection  across  the  x  and  y 
axes  will  extend  the  graphs  to  all  four  quadrants. 

Since  we  can  place  a  conductor  along  any  equipotential  line  we  can 
convert  the  system  to  an  asymmetric  two-plate  transmission  line  with  an 
infinite  ground  plane  at  y  -  0.  As  can  be  seen  from  the  plots  a  low  value  b/a 
results  in  a  minimum  fringing  effect  at  the  edges  of  the'upper  plate. 


5.5  Field  Variation 


The  field  variation  across  an  asymmetrical  tvo-plate  transmission  line 
is  Illustrated  in  fig.  22.  The  values  are  normalized  to  that  of  a 
semi-infinite  system  and  the  dimension  x  is  measured  from  the  edge  of  the 
upper  plate.  As  may  be  seen  the  field  strength  drops  to  50X  at  the  plate 
edge.  Longitudinal  fringing  will  mainly  be  found  at  the  taper  junctions.  For 
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instance,  a  taper  ratio  of  3:1  results  in  a  43;  rise  in  field  strength  at  the 
junction. 


5 . 6  Summary 

If  one  is  prepared  to  tolerate  field  variations  vithin  a  factor  of  two, 
it  is  possible  to  utilize  a  line  twice  the  object  height  with  objects 
separated  by  a  distance  equal  to  their  height  and  work  to  the  edge  of  the 
line. 

For  3031  variations  the  line  height  needs  to  be  3  times  the  object 
height  and  objects  should  be  spaced  apart  by  1%  times  their  height  and  kept  1 
object  height  away  from  the  line  edge. 
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6.0  DESIGN  OF  A  6M  TRANSMISSION  LINE 

6.1  Dimensions 

The  dimensions  of  the  simulator  depend  on  the  size  of  the  objects  in 
the  working  volume.  Typical  large  objects  for  the  Army  are  vehicles,  tanks, 
etc.  and  for  the  Air  Force  small  aircraft  and  helicopters.  The  height  of  these 
large  metal  objects  is  assessed  as  3  and  5m  respectively.  If  field  variations 
of  50Z  can  be  tolerated  the  height  of  the  horizontal  part  of  the  simulator 
should  be  at  least  6m  and  10m  respectively. 

Let  us  start  vith  the  6m  transmission  line.  From  the  graphs  in  fig.  9 
to  20  it  can  be  seen  that  simulators  vith  a  low  b/a  ratio  have  the  best  field 
uniformity  (low  impedance  lines).  Let  us  look  for  a  width  2a  of  6  m,  8  a, 
10  a  and  12  a.  The  ratio  b/a  is  then  2,  1.5,  1.2  and  1.  The  electromagnetic 
field  orientation  and  relative  magnitude  up  to  half  height  and  to  the  edge  of 
the  working  volume  can  be  read  from  the  corresponding  figures.  A  width  of  10m 
appears  to  be  a  reasonable  compromise,  see  the  dashed  line  in  fig.  16.  If 
smaller  objects  are  expected  a  vid'h  of  8  i  is  a  good  choice  as  well  and  it 
reduces  the  apex  angle  4  of  the  input  taper  (shorter  rise  time).  The  useful 
length  of  the  mid-section  is  limited  by  the  reflection  at  the  junction  between 
input  taper  and  mid-section,  see  fig.  23. 

As  mentioned  earlier  the  input  taper  launches  a  TEN  type  wave  vith  a 
spherical  vavefront.  When  this  wave  enters  the  simulator  vorking  volume  a 
reflection  occurs.  This  reflection  in  area  B  adds  to  the  wave  launched  by  the 
taper  sections  at  a  time  corresponding  to  the  clear  time  and  causes  a 
distortion  of  the  wave  onset  and  amplitude.  The  field  in  the  A  portion  of  the 
vorking  volume  shown  in  figure  23,  is  not  distorted  by  the  reflected  vave. 
This  means  that  the  length  of  an  undlstorted  vorking  volume  depends  on  the 
length  of  the  input  taper.  The  input  taper  plays  an  important  role  in  the 
simulator  performance.  A  larger  taper  length  results  in  a  less  spherical 
vavefront,  a  shorter  rise  time,  a  larger  clear  tine  for  reflections  in  the 
direction  of  the  generator  and  in  less  radiation. 
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6.2  Rise  Time 


For  a  length  of  the  working  volume  of  20  m  at  a  working  height  of  3  m  a 
taper  length  of  40  m  is  needed  to  be  outside  the  reflection  area.  The 
spherical  wave  launched  in  the  input  taper  must  be  converted  to  o  plane  wave 
in  the  working  .volume  to  establish  the  waveguide  TEM  propagating  wave.  This 
conversion  takes  place  in  a  time  that  is  equivalent  to  the  propagation  time 
corresponding  to  the  spherical  wave  distance  At,  see  fig.  23: 

At  *  h  (esc  0  -  cot  0  )  (7) 


The  time  associated  with  propagation  over  this  distance  is 
representative  of  the  TEH  wave  rise  time.  Vith  a  taper  length  of  40  m  and  a 
simulator  height  of  6  m,  0  -  8.53°  and  At  -  0.45  m.  The  corresponding  rise 
time  t^  =  At/c  -  1.5  ns.  The  clear  tine  of  the  taper  is  2x40/c  »  267  ns. 


The  total  rise  time  of  the  simulator  +  pulser  is: 

-  h2  *  'P2  <»> 

For  a  total  rise  time  t{  «  5  ns  the  rise  time  of  the  pulser  should  be 

t  i  4.77  ns. 

P 

If  desired  the  input  taper  can  be  used  as  working  volume  with  the 
restrictions  mentioned  earlier.  The  horizontal  vave-guide  acts  then  as  a 
perfect  termination  for  the  taper,  much  better  than  in  the  case  of  the  conical 
plate  simulator.  Horeover  the  radiation  level  is  lover  than  vith  the  CPS. 


The  input  taper  can  be  made  shorter  at  the  sacrifice  of  some  vorking 
height  at  the  end  of  the  vorking  volume.  Vith  a  taper  length  of  30  m  the  rise 
time  of  the  line  t  -  2  ns.  The  clear  time  is  200  ns  and  the  rise  time  for 
the  pulser  should  be  tp  £  4.58  ns.  The  rise  time  measured  in  practical 
situations  at  the  input  of  the  vorking  volume  is  smaller  than  the  calculated 
value.  That  is  because  At  in  fig  23  with  a  vorking  height  of  h/2  is  smaller. 
The  effective  rise  time  of  the  vavefront  should  bo  considered  in  tvo  planes. 
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in  azimuth  as  veil  as  in  elevation.  For  the  calculation  of  the  path 
difference  At  in  azimuth  only  half  of  the  width  of  the  line  should  be  used. 
Vith  a  taper  length  of  30  m  and  a  width  of  10m,  Ai  =  0.41  m  and  t^  =  1.4  ns. 

Since  the  rise  time  of  the  taper  in  both  azimuth  and  elevation  will 
not  contribute  significantly  to  the  overall  rise  time  vithin  the  facility  a 
taper  length  of  30  m  will  be  applied  in  the  further  design. 


6.3  Reflections 


An  abrupt  change  at  the  junctions  between  taper  and  line  will  introduce 
a  capacitive  mis-match  and  result  in  the  mentioned  reflections.  A  radius 
introduced  at  these  points  will  minimize  this  effect  and  improve  the  design. 
The  reflection  at  the  junction  of  the  parallel-plate  section  and  the 
termination  (output)  taper  can  be  reduced  in  tho  same  way  as  with  the  input 
taper.  A  smaller  length  of  the  output  taper  results  in  larger  reflections  and 
shorter  clear  times  of  the  multiple  reflections,  see  fig.  2.  The  parallel- 
plate  section  is  given  the  same  length  as  the  two  tapers  to  make  optimum  use 
of  the  undistorted  area  at  the  ground  plate  (experiments  vith  cables  upon  the 
ground).  This  section  does  not  contribute  to  the  radiation  of  the  line,  see 
appendix  A. 


6.4  Impedance 

The  design  aim  of  a  simulator  transmission  line  is  to  keep  the 
characteristic  impedance  Z as  constant  as  possible  from  the  beginning  to  the 
end  of  the  line  to  prevent  reflections.  This  is  a  k*nd  of  an  art,  since  many 
factors  are  influencing  the  impedance,  some  of  vhich  are  not  very  well 
understood.  The  design  can  be  started  with  calculating  the  approximate  values 
of  the  impedance.  Afterwards  the  impedance  can  be  measured  vith  TDR  methods. 
Vhere  necessary,  corrections  can  be  made  in  the  structure. 

The  derivation  of  the  line  impedance  has  been  extensively  treated  by 
Baum  in  his  Sensor  and  Simulator  Note  SSN  21  (12).  Brovn  et  al  used  this  note 
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to  produce  a  set  of  graphs  and  tables  to  simplify  the  design,  see  tables  1  to 
5  and  figs.  9  to  20  (13).  The  impedance  can  be  defined  as: 

Z  «  fg*Z0  »  (9) 

where  f  is  a  geometric  factor  of  the  particular  line  and  Z  is  the  wave 
impedance  of  free  space.  Factor  f  can  be  read  from  the  tables  1  to  5. 
Parameter  b/a  is  for  asymmetric  two  plate  lines  as  defined  in  fig.  8B,  a  -  1/2 
width  and  b  =  h  (height).  In  our  case  b/a  *  6/5  =  1.2  and  Z  =  0.5‘fg*Zo 

=  98.5  a. 


Fig  21  gives  a  quick  view  of  Z  versus  b/a.  The  generator  and 
termination  dimensions  are  matched  to  those  of  the  line  by  gradual  tapers. 
The  ratio  b/a  applies  for  tht  tapers  as  well. 

It  should  be  noted  from  the  figs.  9  to  20  that  the  field  falls  off 
slowly  with  increasing  distance  from  the  edge  of  the  upper  plate.  Therefore 
the  ground  plane  should  extend  significantly  beyond  the  edges  of  the  upper 
plate  and  intercept  most  of  the  field  for  the  calculations  to  be  valid.  ’When 
the  width  of  the  lover  plate  is  decreased  from  infinity  to  the  same  width  as 
the  upper  plate,  the  impedance  vill  increase.  This  is  shown  in  fig  24.  The 
asymptote  given  at  each  curve  is  the  impedance  for  infinite  width  (14).  In 
practice  the  lover  plate  may  be  in  contact  with  the  earth,  which  is  also  a 
conductor.  This  will  increase  the  effective  width  of  the  lover  plate  and  thus 
decrease  the  line  impedance.  The  effect  of  the  earth  on  the  impedance  is 
negligible  if  factor  a^/a  is  greater  than  three  to  four.  However  that  is  a 
rather  large  ratio.  Therefore  the  earth  is  an  uncertain  factor  in  determining 
the  line  Impedance  in  many  practical  cases.  In  our  case  vith  b/a  *  1.2,  a/h  - 
0.83.  Take  a  ratio  a^/a  •  2.  The  increase  in  line  impedance  is  about  14  Q. 
Dependent  on  the  conductivity  of  the  earth  this  varies  from  14  Q  to  0  Q. 


6.5  Vire  Grids 


Up  to  this  point,  solid  conducting  plates  were  considered  to  form  a 
transmission  line.  Por  large  transmission  lines,  like  this  one,  the  solid 
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plates  should  be  replaced  with  wires  installed  parallel  to  the  direction  of 
propagation  to  avoid  mechanical  or  structural  problems  (weight,  wind,  snow, 
rain,  etc.). 

Wire  grius  will  distort  the  field  near  the  vires  and  will  increase  the 
line  impedance.  The  field  in  the  vicinity  of  a  wire  is  illustrated  in  fig.  25 
(12).  The  leakage  of  the  field  through  the  grid  has  the  effect  of  shifting 
the  grid  to  give  an  increase  in  the  effective  plate  spacing.  This  will  have  a 
corresponding  eftect  on  the  line  impedance.  This  shift  of  the  grid  is  given 
by: 

»  ■  n  l"  te)  (10) 

where  s  is  the  wire  spacing  and  r  is  the  wire  radius. 

If  the  transmission  line  is  removed  from  the  effects  of  the  conducting 
earth,  a  Ah  factor  will  also  be  required  when  the  lower  plate  consists  of 
parallel  vires. 

The  wire  spacing  must  be  small  compared  with  the  wavelength  of  the 
highest  frequency  to  be  propagated  down  the  line.  Vith  wider  wire  spacing, 
the  high  frequencies  will  radiate  energy  between  the  vires.  This  radiation 
will  cause  a  loss  of  high  frequency  content  in  the  wave  and  will  maintain  a 
spherical  wavefront  for  high  frequencies  within  the  working  volume  (15). 

V<th  a  wire  spacing  s  -  40  cm  and  wire  radius  r  -  4  mm,  Ah  -  17.6  cm 
for  the  upper  plate,  and  0  to  17.6  cm  for  the  lover  plate,  dependent  on  the 
conductivity  of  the  earth.  Assume  the  effect  of  the  earth  is  negligible,  then 
the  effective  plate  spacing  is  h  ♦  2*Ah  -  6  ♦  0.352  -  6.352  a.  Factor  b/a  - 
6.352/5  -  1.27,  and  f  -  0.5383,  Z  -  0.5-f  -  101.5  Q.  With  a  lower  plate 
width  two  times  the  upper  plate  width  14  Q  has  to  be  added:  Z  -  101.5  +  14  > 
115.5  Q.  The  wire  spacing  s  -  40  cm  corresponds  with  a  frequency  c/0.4  - 
750  MHz  which  is  high  compared  to  the  highest  frequency  to  be  propagated,  see 
fig.  4. 


The  accuracy  of  Impedance  suggested  in  the  tables  is  partly  cancelled 


i 
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by  the  unknown  effect  of  the  earth  upon  the  impedance.  This  uncertainty  can 
be  decreased  by  connecting  the  perimeter  of  the  lower  plate  to  ground  rods. 
The  impedance  of  the  line  can  be  measured  with  a  reflectometer ,  an  instrument 
which  measures  the  pulse  impedance  by  reflection  techniques.  It  shows  the 
level  and  the  location  of  the  reflection. 


6.6  Field  Enhancement 


The  replacement  of  solid  plates  by  a  grid  of  wires  has  a  second 
consequence;  a  field  enhancement  around  the  wires,  see  fig.  25.  Shannon  has 
derived  an  expression  for  the  field  enhancement  (16).  For  a  grid  of  parallel 
wires  with  the  field  applied  to  one  side  the  field  enhancement  factor  is: 


f  -  1/2  +  1/12-  pJ!)  + 


(ID 


for  d/s  «  0.1,  where  d  is  the  wire  diameter  and  s  the  wire  spacing.  Fig.  26 
shows  f  versus  d/s. 


Another  field  enhancement  exists  at  the  edge  of  the  upper  plate  or  wire 
grid.  In  our  case,  with  b/a  »  1.3,  see  fig.  17.  For  the  geometry  of  a  semi¬ 
infinite  plate  with  a  radiused  edge  Shannon  gives  as  field  enhancement  factor: 

f  -  for  r/h  «  1  (12) 
n  r 


where  r  is  the  edge  radius  and  h  the  height  of  the  line.  Fig.  27  shows  f 
versus  r/h.  lr.  both  cases  the  local  field  is  clearly  well  in  excess  of  the 
general  field  and  may  in  fact  exceed  the  breakdown  strength  of  air.  A 
discharge  in  the  air,  or  a  spark  to  a  conductor  may  be  the  result. 


6.6.1  Breakdown  of  the  Air 


The  breakdown  strength  of  the  atmosphere  is  generally  given  as 
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25  kV/cm.  The  DC  breakdown  strength  of  a  spark  gap  with  air  in  a  uniform 
field  is  given  as: 


E 


24. 6p  + 


6.7  5 

u 


kV/cm 


(13) 


where  p  is  the  pressure  in  atmospheres  and  d  the  gap  spacing  in  cm.  With 
p  =  1  and  d  =  lf  E  =  31.3  kV/cm. 

For  short  pulses  the  breakdown  field  strength  can  be  larger,  as  a 
result  of  the  time-lag  to  initiate  a  Townsend  avalanche  and  to  form  a  spark 
channel.  It  takes  time  to  form  a  streamer.  For  instance,  with  a  streamer 
velocity  v  *  l(r  m/s  a  pulse  length  of  at  least  10  ys  is  needed  to  cause  a 
spark-over  at  a  distance  of  I  m  (18). 

Let  us  look  to  the  design  of  our  grid  of  wires.  It  has  advantage  to 
construct  this  grid  as  light  and  as  strong  as  possible,  for  instance  apply  not 
too  many,  small  diameter,  stranded  wires.  With  a  diameter  d  ■  3  mm  and 
spacing  s  =  1  m,  d/s  =  3*10~^.  From  fig.  26  it  can  be  read  fhat  the  field 
enhancement  factor  £  =  110.  With  a  main  field  of  0.5  kV/cm  the  local  field  at 
the  wires  is  110  x  0.5  *  55  kV/cm.  That  is  too  large.  Assuming  a  breakdown 
field  strength  of  air  of  25  kV/cm  and  applying  a  safety  factor  of  2,  the 
maximum  allowable  field  enhancement  factor  is: 

.  12.5  kV/cm 

ra  "  U.TkWd  * 

_2 

This  value  corresponds  in  fig.  26  with  a  ratio  d/s  -  1.3*10  .  When 

choosing  a  wire  spacing  s  -  40  cm,  to  prevent  radiation,  it  follows  that 

d  »  5.2  mm.  It  is  recommended  to  choose  a  field  strength  in  the  facility  that 

is  a  little  larger  than  the  required  50  or  65  kV/:.i,  say  Ej  >  75  kV/m.  With  an 

E,  >  75  kV/m,  f  »  17  and  the  diameter  d  ■  8  mm. 
i  m 

The  fie'.d  enhancement  at  the  edge  of  the  grid  of  wires  can  be  approxi¬ 
mated  with  fig  27  (valid  for  solid  plates).  With  a  wire  radius  r  -  4  mm  the 

-4 

ratio  r/h  -  6.7*10  ,  and  f  ■  22.  lh»  enhanced  field  at  the  edge  is  22  x  0.75 

>  16.5  kV/cm.  This  means  the  wire  radius  should  ot  be  less  than  4  mm. 
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In  the  taper,  the  height  h  is  decreasing,  causing  a  lower  field 
enhancement  factor  f  and  a  larger  field  strength  E^  in  the  taper.  Table  6 
shows  the  edge  enhanced  field  Eg  of  the  taper  with  a  facility  field  strength 
Ej  =  75  kV/m  or  0.75  kV/cm  and  a  wire  radius  r  =  4  mm.  It  appears  that  Ee  is 
rather  large,  particularly  in  the  lower  height  regions. 


6.6.2  Larger  Edge  Radius 

To  be  sure  that  no  breakdown  of  the  air  will  occur  the  edge  radius 

should  be  increased.  With  a  factor  f  =  3  the  field  strength  E  is  13.5  kV/m, 

-2  6 
the  ratio  r/h  =  4.5-10  and  r  =  45  mm.  This  can  be  simulated  by  some  thinner 

wires  placed  at  the  circumference  of  a  half  circle  with  radius  r  =  4.5  cm,  see 

fig.  28. 


An  advantage  of  this  construction  is  that  later  on,  the  transmission 
line  can  also  be  used  for  er.do-atmospheric  simulations.  A  disadvantage  is 
that  the  transmission  line  is  less  flexible.  However,  if  a  reinforced 
concrete  pad  is  needed  to  allow  heavy  objects  upon  the  ground  plate  the 
facility  is  more  or  less  stationary  as  it  is.  In  this  case  aluminum  wires 
with  a  diameter  of  8  mm  are  recommended. 


6 . 7  Detailed  Design 

Now  that  the  dimensions  and  the  construction  of  the  facility  regarding 
the  given  boundary  conditions  are  outlined,  a  more  detailed  design  of  the 
transmission  line  can  be  started. 


6.7.1  Horizontal  Section 


The  separation  of  the  upper  and  lower  grid  of  wires  is  6  m.  The  upper 
grid  has  a  width  of  10  m  and  the  lower  20  m,  see  figs.  29  and  30.  They  are 
constructed  of  aluminum  wires  with  diameter  d  -  8  mm  passed  through  plywood 
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spreaders  which  stiffen  and  shape  the  line.  The  spacing  between  the  wires  is 
40  cm.  The  upper  grid  has  26  wires  and  the  lower  51  wires.  At  the  edges  of 
the  upper  grid  5  extra  wires  with  diameter  d  =  3  mm  are  forming  a  rounded  edge 
with  a  radius  r  =  5  cm.  The  lower  grid  is  at  the  periphery  connected  to  small 
ground  rods  to  minimize  variations  in  the  influence  of  the  ground  and  to 
decrease  the  impedance. 

The  effect  of  the  rounded  edges  upon  the  line  impedance  has  to  be 
compensated.  Suppose  the  rounded  edges  can  be  replaced  by  flat  plates  with  a 
width  of  approximately  l/4nd  ~  8  cm.  To  keep  the  same  impedance  the  flat  part 
of  the  upper  grid  should  have  a  width  of  10  m  -  (2  x  8  cm).  The  wire  spacing 
s  »  40  cm  is  then  a  little  smaller. 


6.7.2  Input  Taper 

It  is  assumed  that  the  height  of  the  generator  output  is  about  1  m  and 
the  peak  output  voltage  is  450  kV.  Table  7  shows  the  enhanced  field  at  the 
edges  and  between  the  wires  at  several  heights  of  the  taper.  The  maximum 
values  at  the  generator  connection  are  13  kV/cm  and  17  kV/'cm  respectively. 
The  field  strength  between  the  wires  can  be  reduced  by  constructing  the  first 
part  of  the  taper  of  solid  plate.  The  connection  to  the  generator  should  be 
given  round  edges  in  all  directions  to  prevent  additional  field  enhancements. 
The  remaining  maximum  edge  field  strength  is  then  13  kV/cm,  a  factor  2  below 
the  atmospheric  breakdown  level.  This  safety  factor  is  a  minimum,  since  an 
incorrect  adjustment  of  the  generator  can  double  the  output  voltage,  see 
section  6.8.2. 

With  these  precautions  complex  structures  with  gas  tight  compartments 
filled  with  sulphur  hexafluoride  gas  (SF^)  can  be  avoided.  To  close  the 
generator  housing  and  to  electrically  isolate  the  taper  from  the  housing  a 
perspex  interface  can  be  used  with  a  thickness  of  about  15  mm.  The  minimum 
distance  from  conductive  parts  of  the  construction  should  be  at  least  1  m. 

In  table  8  the  impedance  of  the  taper  at  several  places  is  given,  see 
fig.  31.  In  the  calculations  an  effective  height  is  assumed  of  b'  -  h  ♦  2 Ah. 
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If  a  part  of  the  upper  or  lower  grid  is  made  of  solid  plate  the  correction 
with  Ah  or  2 Ah  is  not  needed.  The  impedance  of  the  taper  is  decreasing  a 
little  with  decreasing  height.  This  can  be  corrected  by  increasing  the 
height.  The  correction  in  fig.  31  is  given  in  centimetres. 

In  practice  the  impedance  is  a  few  ohms  larger  because  of  the  limited 
width  of  the  groundplate.  The  interface  between  the  solid  part  and  the  wire 
part  of  the  taper  should  be  made  gradual,  e.g.  with  wire-mesh.  At  the  top  of 
the  taper  the  reflections  can  be  reduced  by  using  a  small  elevation  angle 
0  =  11.31°  and  by  rounding  of  the  shape  of  the  junction. 


6.7.3  Output  Taper 

* 

The  output  taper  has  the  same  dimensions  and  construction  as  the  input 
taper.  The  effective  rise  time  is  2  ns  in  elevation  and  1.4  ns  in  azimuth. 
The  main  task  of  the  output  taper  is  gradually  connecting  the  horizontal 
parallel-plate  section  to  its  load.  The  transmission  lihe  is  distributed 
loaded  with  a  construction  that  forms  an  integral  part  with  the  taper.  The 
dimensions  of  the  connection  to  the  load  are  determined  by  the  terminating 
construction. 

6.7.4  Termination 


The  aim  of  a  good  design  is  to  prevent  reflection  of  the  guided  wave 
back  into  the  working  volume.  The  output  taper  guides  the  wave  to  a 
relatively  small  resistive  load,  whose  impedance  is  equal  to  the 
characteristic  impedance  of  the  transmission  line  structure  and  is  "invisible" 
for  the  arriving  wave  front.  The  best  way  to  prevent  reflections  is  to  make 
the  transmission  line  infinitly  long.  This  can  be  simulated  by  distributing 
the  load  over  a  finite  length  and  make  it  a  part  of  the  taper. 
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Other  design  aims  are: 

-  the  termination  should  be  resistive  for  all  frequencies  of  interest, 
see  Fig  4 

-  the  termination  should  be  able  to  handle  a  peak  voltage  of  450  kV 

-  the  resistance  should  be  adjustable  between  100  and  110  Q 

A  good  choice  for  the  resistive  material  is  the  carbon  type  resistor. 

. » 

It  consists  of  a  ceramic  tube  with  a  thin  carbon  layer  on  the  outside.  A 
suitable  choice  xs  Morganite,  type  701  with  a  low  inductance  and  ar.ti-track 
bands  to  avoid  spark-overs.  The  length  is  250  ram,  the  diameter  25  mm. 


With  the  anti-track  bands  the  maximum  allowable  voltage  per  resistor  is 
30  kV.  To  handle  450  kV,  15  resistors  in  series  are  needed.  With  20 
resistors  in  series  a  safety  factor  of  1.3  is  created  and  the  total  length  is 
in  agreement  with  the  minimum  height  of  1  m,  see  Pig.  32.  The  inductance  of 
20  resistors  in  series  can  be  reduced  and  the  termination  can  be  distributed 
over  the  whole  width  by  connecting  a  number  of  resistor-strings  in  parallel. 


With  16  resistor-strings  in  parallel  and  82  fi  per  resister,  the  total 
value  is: 


20  x  82 

t  “  T5 - 


102.5  a 


By  disconnecting  one  of  the  strings  the  total  value  ■  109.3  a. 

Damage  by  over-voltage  across  one  of  the  resistors  can  be  prevented  by 
connecting  the  16  strings  not  only  at  the  two  ends  but  also  at  each  resistor 
station;  thus  forming  equipotential  lines,  see  fig.  33.  By  removing  one  or 
more  resistors  fine  adjustments  can  be  made  (don't  overstress  the  other 
resistors). 

Since  the  resistive  network  forms  a  part  of  the  taper  and  has  a  space 
dependent  resistive  value,  the  impedance  of  that  part  of  the  taper  (formed  by 
the  terminating  network)  should  be  matched  to  the  resistive  value.  This  can 
be  done  by  increasing  the  width  of  the  network  to  infinity  at  zero  height.  In 
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practice  a  width  two  times  that  at  1  m  height  is  sufficient.  The  ground  plate 
is  made  of  solid  material.  The  rounded  edges  of  the  taper  continue  in  the 
terminating  network  by  two  resistor  strings  at  each  side.  Fig.  34  shows  some 
construction  details  of  the  edges  and  in  figs.  35  and  36  a  suggestion  for  the 
wooden  frame  work  is  given.  A  detailed  field  enhancement  plot  of  the  edges  can 
be  made  with  the  "resistive  paper  and  silver  paint"  method.  With  the  proposed 
dimensions  of  the  terminating  network  the  enhancement  is  up  to  the  breakdown 
limi t . 

The  peak  voltage  across  the  resistive  network  V  =  450  kV.  With  a  value 
R  ~  100  2  the  peak  current  i  ^  =  4.5  kA.  For  a  half  width  time  t^  =  0.7t  = 
200  ns,  C  =  2.86  nF  and  W  =  1/2 ’CV^  =  290  J.  The  t^  of  an  endo-atmospheric 
EMP  is  about  190  ps  (exponential  decay)  and  W  =  27.5  kJ.  In  both  cases  the 
dissipated  energy  is  far  below  the  maximum  rating. 


6.8  Pulse  Generator 

The  pulse  generator  (also  called  pulser)  loaded  with  about  100  Q, 
should  have  a  peak  output  voltage  of  450  kV  and  a  time  history  similar  to  that 
of  the  exo-atmospheric  EMP,  see  fig.  3.  The  field  strength  E(t)  has  the  same 
shape  versus  time  as  the  output  voltage  of  the  pulser.  The  peak  value  is 
450  kV/6m  *  75  kV/ia. 

High  voltage  generators  have  relatively  large  dimensions,  even  when 
they  generate  this  voltage  for  a  short  time.  As  a  consequence  they  contain 
too  much  inductance  to  deliver  high  output  voltages  with  a  rise  time  in  the  ns 
range.  High  voltages  together  with  short  rise  times  car  be  obtained  by  a 
combination  of  a  Marx  generator  and  a  pulse  forming  network  like  a  peaking 
capacitor  or  a  transfer  capacitor.  In  this  case  the  Marx  generator  should 
have  an  output  voltage  of  about  450  kV  +  20%  -  540  kV.  Ask  for  a  600  kV  Marx 
to  prevent  an  unpleasant  surprise  during  the  acceptance  test. 


-  26  - 


The  basic  peaking  capacitor  circuit  is  schematically  given  in  fig.  37. 

C  is  ihe  total  capacitance  of  the  Marx  generator  and  L  is  its  inductance 

n  m 

lumped  together.  Cp  is  the  peaking  capacitance  and  is  the  load  resistance. 

Initially  C  is  charged  to  V  volts,  the  Marx  voltage.  After  switch  S. 

no  i 

is  closed  the  peaking  capacitor  will  be  charged.  The  value  of  is  chosen 
so  that  when  the  voltage  across  it  reaches  VQ,  the  current  through  it  is 
V  /R  .  At  this  time  Sn  is  closed  and  the  voltage  across  R  and  the  current 

Ol  4  t 

through  R  become  immediately  VQ  and  Vq/Ri  respectively.  In  practice  the  rise 
time  of  the  voltage  and  current  are  greater  than  zero  because  C  ,  Rt  and  the 
interconnecting  conductors  contain  some  inductance. 


The  peaking  capacitor  is  a  relatively  small,  high-voltage,  low- 
inductance  capacitor.  The  capacitance  is  large  enough  to  build  up  the  current 
from  the  Marx  generator  to  a  level  required  by  the  load.  When  the  generator 
current  has  build  up,  the  load  is  connected  and  the  current  is  transferred  to 
the  load,  see  fig.  38. 


(14) 

(15) 


i,(t)  is  the  current  through  L  and  V  ft)  is  the  voltage  across  C 

u  m  p  p 


At  t 
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V 

This  value  should  be  equal  to 

i 


This  means  R  =  u  L  =  JL  /C 
i  o  m  m 


or  C  =  —x 
P  n  2 


(16) 


The  decay  of  the  pulse,  after  Sg  is  closed,  is  determined  by  the 
capacitance  (C^  +  Cp)  and  R,.  The  required  half  width  time  0.7t  *  200  ns, 

,-9 


(C  +  C  )  =  £00^0  =  2  g6  nF> 

v  m  p'  0.7  •  P. 


For  C  =  0.1-C  ,  C  ~  290  pF  and  L  »  2.9  yH. 
p  in  p  n 

The  calculation  of  Cp  is  an  approximation,  since  some  stray  inductance 
and  parasitic  capacitance  is  neglected.  The  peaking  capacitor  can  be  made 
fairly  small,  with  a  low  inductance.  With  a  good  design  of  Cp  the  rise  time 
of  the  pulser  (Marx  generator  +  peaking  circuit)  is  mainly  defined  by  the 
series  inductance  of  the  spark-gap  S2» 


6.8.2  Spark  Gap 

Spark  gap  S^  in  fig.  37  is  the  spark  gap  in  one  of  the  stages  of  the 
Marx  generator  when  they  are  switched  into  series.  It  is  a  low  voltage 
(30  kV)  switch  that  does  not  play  an  important  role  in  the  total  circuit.  The 
main  spark  gap  S2  should  ideally  have  no  inductance,  however  tne  rise  time  is 
controlled  by  two  factorss  the  inductance  of  the  spark  channel  and  the  rate 
of  current  growth  during  the  resistive  phase,  so  called  because  of  the  energy 
used  to  heat  and  expand  the  plasma  channel.  The  inductance  is  related  to  the 
gap  distance  and  the  resistive  phase  to  the  density  of  the  applied  gas. 

The  gap  distance  can  be  reduced  by  applying  a  gas  with  a  larger 
electric  strength  than  air,  e.g.  SFg  with  a  factor  2  to  2.3  related  to  air  and 
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by  increasing  the  gas  pressure.  To  further  decrease  the  inductance  of  the 
spark  gap  the  dimensions  should  be  made  as  small  as  possible. 

In  our  case  the  voltage  across  the  spark  gap  is  450  kV.  To  prevent 
sparkovers  upon  the  outside  of  the  spark  gap  housing,  field  enhancements 
should  be  avoided  and  the  spark  gap  can  be  placed  in  oil,  together  with  the 
peaking  capacitor. 

The  time  the  main  spark  gap  fires  determines  the  shape  of  the  output 
voltage  of  the  pulser,  see  fig.  39  (20).  In  extreme  cases  the  output  voltage 
can  be  two  times  the  desired  value  (2-Vo  in  fig.  38).  The  firing  of  the  spark 
gap  can  be  triggered.  The  time  for  self  firing  can  be  adjusted  with  the  gap 
distance  and  the  gas  pressure. 


6.8.3  Voltage  Divider 


Vhen  the  peaking  capacitor  and  the  spark  gap  vith  trigger  capacitor  are 
placed  in  oil  to  keep  the  dimensions  small,  also  a  coaxial  voltage  divider  can 
be  combined  vith  it  to  monitor  the  input  voltage  of  the  transmission  line. 


6.9  Electrical  Specifications  of  the  Pulser 


Vhen  ordering  a  pulser  the  following  specifications  are  important! 

-  maximum  output  peak  voltage:  500  kV 

-  minimum  output  peak  voltage:  100  kV 

-  rise  time:  5  ns  or  less 

-  load:  resistive,  100-110  Q 

-  tall  time:  e-fold  time  286  ns 

-  prepulse:  <  IX  of  peak 

-  oscillations:  <  52  of  momentary  amplitude 

-  pulse-to-pulse  repeatability:  £  22  deviation 

-  output  polarity:  positive  and  negative 

-  charging  time:  1  PPM  at  500  kV 

-  lifetime  1 

-  spare  parts  >  (19) 

-  definitions  ) 
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In  this  report  only  the  electrical  properties  of  the  simulator  have 
been  considered.  Hechanical  specifications  such  as  the  matching  of  the  pulser 
output  to  the  transmission  line  are  also  very  important. 
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7.0  DESIGN  OF  A  10m  TRANSMISSION  LINE 

If  larger  objects  have  to  be  investigated  or  the  simulator  should  not 
have  a  fixed  base,  it  is  worth  while  to  look  at  which  of  the  properties  can  be 
sacrificed  for  larger  dimensions  and  more  flexibility. 


7.1  Dimensions 


Let  us  use  a  helicopter  as  an  example.  The  largest  height  is  assessed 
as  5.2  m  and  the  length  as  19  m.  For  exposing  the  helicopter  in  all  direc- 
tions  a  working  volume  is  needed  of  about  i  x  v  x  h  •  20  x  20  x  10  m  .  If 
field  variations  of  50Z  are  acceptable  the  simulator  height  can  be  limited  to 
10  m  and  the  width  to  20  m.  Vhen  the  same  pulser  of  the  6m-  line  is  applied 
with  an  output  peak  voltage  of  450  kV  the  peak  field  strength  is  450  kV/lOm  - 
45  kV/m.  Vith  a  height  of  10  ■  and  width  of  20  m,  the  ratio  of  b/a  -  1.  The 
uniformity  and  orientation  of  the  field  can  be  read  from  fig.  14.  They  are 
better  than  from  the  6  m-line. 


7.2  Rise  Time 


To  have  the  same  performance  as  the  6  a  line  the  input  taper  should 
have  a  length  of  5h  that  means  5  x  10  -  50  a.  That  is  too  long  to  be  flexible 
and  there  is  not  always  place  for  a  simulator  vith  large  tapers.  Therefore 
some  distortion  caused  by  reflections  should  be  accepted. 

Let  us  try  a  length  of  20  or  30  a.  In  this  case  the  rise  time  in 
horizontal  and  vertical  direction  is  the  same.  Vith  20  a  length  tr  -  13  ns 
and  with  30  m  tr  -  5.4  ns.  Since  the  responses  of  exposed  equipment  are 
directly  proportional  to  the  rise  time  (the  first  derivative  of  the  field)  no 
further  concession  rhould  be  Bade  in  this  direction.  Ve  Bake  a  choice  for  the 
30  a  taper.  Another  advantage  is  a  larger  undlstorted  working  volume. 
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7.3  Impedance 

The  impedance  of  the  line  with  infinite  width  of  the  ground  plate  is 
89.1  Q.  Some  of  the  fine  uniformity  of  the  field  in  fig.  14  has  to  be 
sacrificed  for  more  flexibility.  Ve  give  the  ground  plate  the  same  width  as 
the  upper  plate.  With  a  ratio  a^/a  =1  in  fig.  24  the  impedance  will 
increase  to  about  115  2.  The  impedance  can  be  made  smaller  by  making  the 
ground  mo  a  conductive  and/or  to  connect  the  ground  plate  to  the  ground  via 
ground  rods. 

Instead  of  changing  the  resistive  value  of  the  termination,  the 
impedance  of  the  line  can  be  adjusted  by  changing  the  width  or  the  height  of 
the  upper  plate,  e.g.  with  a  height  of  9  m  the  impedance  is  115  2  -  (89.1  - 
83.9)  «  109.8  2  and  the  field  strength  is  50  kV/ra.  The  line  should  be 
flexible! 


7.4  Wire  Grids 

The  replacement  of  plates  by  wires  has  the  additional  advantage  that 
the  height  and  width  of  the  upper  grid  csn  be  changed  easier.  Ve  choose 
stranded  steel  wires  with  a  diameter  d  »  3  mm  and  a  spacing  s  *  1  m.  A 
wavelength  A  -  1  ■  corresponds  with  a  frequency  of  f  -  c/A  -  300  HHz.  Some 
radiation  of  high  frequencies  can  be  expected,  with  a  little  more  spherical 
wavefront  compared  to  the  6  m  line.  Since  also  the  rise  time  is  increased 
from  2  ns  to  5.4  ns  it  is  not  too  bad. 

The  effective  height  of  the  line  is  increased  by  the  wires  with 
Ah  -  74  cm  for  the  upper  grid  and  74  cm  for  the  lower  grid.  This  2* Ah  can  be 
limited  to  74  cm  by  applying  chicken  wire  at  the  ground.  This  has  the 
additional  advantage  that  it  is  more  flexible  and  can  be  extended  to  a  larger 
area  (lower  impedance). 
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7.5  Field  Enhancement 


7.5.1  Horizontal  Section 


Between  the  vires:  d/s  -  3*10  and  f  -  110.  With  a  facility  field 
Ej  ■  0.5  kV/cm  the  local  field  is  110  x  0.5  *  55  kV/cm.  That  is  too  large! 
The  breakdown  voltage  of  the  atmosphere  is  about  25  kV/cm.  Let  us  try  a 
spacing  s  *  50  cm,  d/s  -  B’lO-^,  f  «  55  and  Ev  -  55  x  0.45  *  25  kV/cm. 

At  the  edges:  r/h  -  1.5* 10“^,  f  *  46,  Ee  -  23  kV/cm.  In  this  case  the 
enhanced  field  between  the  vires  and  at  the  edges  is  about  the  same  as  the 
breakdown  voltage  of  the  air.  For  a  breakdown  a  minimum  pulse  length  is 
needed.  If  we  use  the  half  width  time  of  the  EHP,  t^  -  200  ns  and  assume  a 
streamer  velocity  v  ■  10  cm/us  a  breakdown  can  be  expected  at  a  distance  of 
2  cm.  Therefore  a  minimum  distance  of  at  least  5  cm  should  be  used.  Also 
corona  should  be  avoided.  Corona  is  dependent  on  factors  such  as  humidity, 
air  pressure,  polarity,  etc.  If  only  the  exo-atmospheric  EHP  is  simulated 
with  a  short  pulse  length  and  any  other  field  enhancement  can  be  avoided,  it 
may  be  possible  to  go  to  the  limits.  Another  possibility  is  applying  a  lover 
output  voltage  of  the  pulser. 


7.5.2  The  Input  Taper 

Assume  that  the  height  of  the  pulser  output  is  1  m.  Let  us  laok  at  the 
field  enhancement  at  this  point  because  the  normal  local  field  strength  is 
here  450  kV/m  -  4.5  kV/cm. 

_2 

Between  the  vires:  the  spacing  s  -  5  cm,  d/s  -  6*10  ,  f  -  6.3  and  the 

local  field  Ey  -  6.3  x  4.5  -  28.4  kV/cm. 

At  the  edges:  r/h  -  1.5-10'^,  f  •  14.5  and  B#  -  14.5  x  4.5  - 
65.3  kV/cm. 


Table  9  shows  the  enhanced  field  between  the  vires  and  at  the  edges  for 
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a  facility  field  strength  =  50  kV/m.  The  wire  spacing  s  =  50  cm  and  the 
wire  diameter  d  =  3  mm.  From  a  taper  height  of  4  m  down  to  1  m  the  local 
field  Eg  increases  from  36  to  72  kV/cm.  This  can  be  reduced  by  increasing  the 
radius  of  the  edges,  for  instance  with  thin  wire  mesh  (chicken  wire).  A 
radius  of  1  cm  is  enough,  gradually  decreasing  with  height.  The  same  wire 
mesh  can  be  used  to  short  circuit  the  wires  in  the  lower  part  of  the  taper, 
thus  reducing  Ev»  I*,  is  also  possible  to  use  thicker  wires  at  the  edges,  also 
for  mechanical  reasons.  With  no  further  field  enhancements  the  local  field  is 
everywhere  around  25  kV/cm.  And  that  is  the  limit! 


7.6  Detailed  Design 


7.6.1  Horizontal  Section 


The  upper  grid  has  a  width  of  20  m  at  a  height  of  10  m,  The  wires  are 
spaced  at  50  cm  and  have  a  diameter  of  3  mm.  The  grid  has  2C/0.5  *  1  »  41 
wires.  The  length  is  dependant  on  the  dimensions  of  the  object  under  test  and 
the  length  of  the  input  taper.  In  fig.  40,  a  maximum  length  of  30  m  has  been 
chosen.  The  lower  plate  can  be  made  of  wire-mesh,  connected  to  earth  at  the 
periphery  for  lower  impedance  and  better  field  uniformity. 


7.6.2  Input  Taper 

The  field  enhancement  shown  in  table  9  can  be  reduced  by  applying  solid 
metal  plate  or  wire-mesh  at  the  taper  end,  near  the  pulser.  To  reduce  E@  at 
1  m  height  from  72.5  to  30  kV/cm,  the  field  enhancement  factor  f  should  be 
made  smaller,  e.g.  30/72.5  x  14.5  -  6.  In  fig.  27,  a  ratio  r/h  *  9*10  can 
be  read,  which  results  in  a  radius  of  9  mm.  The  effective  width  is  increased 
by  this  radius  with  2'6a,  approximated  by  l/4nd  ■  1.57  cm  for  each  edge. 

Table  10  shows  the  fine  adjustment  for  a  constant  impedance  Zx  -  91  S. 
The  taper  can  be  lifted  up  a  few  centimeters  at  the  pulser  end. 
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The  exact  value  of  Z^  is  not  so  important,  only  the  constant  value. 
When  the  lover  plate  of  wire  mesh  is  in  contact  with  a  good  conductive  ground 
the  impedance  Z^  of  the  line  is  estimated  to  be  91  +  13  =  104  S. 


7.6.3  Output  Taper 

A  width  of  20  m  is  relatively  large  for  one  taper.  Therefore  two 
tapers  are  designed  to  terminate  the  line.  Advantages  are:  a  shorter  rise 
time  in  azimuth,  less  fringing  at  the  junctions,  better  uniformity  of  the 
field  and  a  simpler  construction  of  the  terminating  network. 


7 .  (y  4  Termination 

The  minimum  height  for  450  kV  should  be  1  m.  The  length  of  the 
terminating  networks  is  then  about  1/10  x  30  »  3  m.  The  network  is  forming  a 
part  of  the  taper.  The  width  is  1/10  x  10  »  1  m  at  the  taper-network 
Junction,  see  fig.  42,  To  match  the  impedance  of  the  line  to  the  local 
resistive  value  of  the  network,  the  width  is  increasing  towards  the  ground 
connection.  The  connection  to  ground  should  have  a  low  impedance.  In  fig.  41 
the  resistors  are  connected  via  a  metal  profile  to  the  reinforcement  of  the 
concrete  block.  The  resistors  are  confined  in  fibre  glass  reinforced  tubes. 
The  impedance  of  the  line  Z(  •  104  fi.  Each  network  should  therefore  have  a 
value  of  208  fi,  or  3  x  208  -  624  fi  per  resistor  string.  Using  the  same  type 
of  resistors  as  in  the  6  m  line  a  total  of  12  resistors  of  52  fi  in  series  are 
needed.  Vith  a  maximum  voltage  per  resistor  of  30  kV,  the  total  voltage  that 
can  be  handled  is  360  kV.  Vith  a  short  pulse  length  a  voltage  overstress  vith 
25*  to  450  kV  can  be  tried.  Otherwise  the  fibre  glass  tubes  can  be  filled 
with  SP^. 
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8.0  CONCLUSIONS 

The  choice  between  a  parallel  plate  and  a  conical  plate  simulator  has 
been  made  in  favour  of  the  first  one.  The  advantages  of  both  are  more  or  less 
combined  by  applying  a  long  input  taper. 

Two  designs  for  a  bounded  wave  simulator  have  keen  made,  one  with  a 
fixed  6  m,  and  another  with  a  more  flexible,  10  m  transmission  line.  For  the 
first  one  a  conservative  design  has  been  used.  The  larger  working  volume  and 
the  better  flexibility  of  the  second,  more  relaxed  design  have  been  obtained 
at  the  cost  of  a  larger  rise  time,  lower  field  strength  and  a  higher  radiation 
level  (see  table  11  and  Appendix  A).  With  the  procedure  given  in  this  report 
any  design  in  between  these  two  can  be  made. 

*- 

No  efforts  have  been  made  for  the  design  of  a  complete  simulator 
(pulser  and  transmission  line)  since  that  was  out  of  the  scope  of  this  task. 
However  the  electrical  and  mechanical  matching  of  a  future  pulser  to  the 
transmission  line  can  be  made  in  a  flexible  way. 

ACKNOWLEDGEMENTS 

I  like  to  thank  Dr.  Bev  G.  Young  and  Dr.  Satish  Kashyap  for  their 
encouragement  and  h.lp  in  preparing  this  report  and  M.A.  O'Keefe  for  typing 
the  manuscript. 


-  36  - 


9.0  REFERENCE, 

1.  NATO  Documents  on  Nuclear  Survivability  Criteria  for  Armed  Forces 
Material  and  Installations. 

2.  The  Electromagnetic  Field  in  an  EMP  Simulator  at  a  High  Frequency, 
R.W.P.  King,  D.J.  Blejer.  IEEE-EMC-21,  no. 3,  Aug  1979. 

3.  Standing  Waves  and  Notches  in  an  EMP  Simulator  and  Their  Reduction, 
R.W.P.  King,  D.J.  Blejer,  T.T.  Wu,  IEEE,  EMC-23,  no.  2,  Hay  1981. 

4.  The  Rhombic  EMP  Simulator,  H.M.  Shen,  R*.W.P.  King,  IEEE,  EHC-24,  no.  2, 
May  1982. 

5.  Experimental  Investigation  of  the  Rhombic  EMP  Simulator:  Comparison 
with  Theory  and  Parallel-Plate  Simulator,  H.M.  Shen,  R.W.P.  King,  IEEE, 
EMC-24,  no.  3,  Aug  1982. 

6.  An  experimental  Investigation  of  the  Rhombic  EMP  Simulator  Under  Pulse 
Excitation,  H.M.  Shen,  R.W.P.  King,  T.T.  Wu,  IEEE,  EMC-25,  no.  1, 

Feb  1983. 

7.  Theoretical  Analysis  of  the  Rhombic  Simulator  Under  Pulse  Excitation, 
H.M.  Shen,  R.W.P.  King,  T.T.  Wu,  IEEE,  EMC-25,  no.  1,  Feb  1983. 

8.  An  Experimental  Investigation  of  the  Parallel-Plate  EMP  Simulator  with 
Single-Pulse  Excitation,  H.M.  Shen,  R.W.P.  King,  T.T.  Wu,  IEEE,  EMC-25, 
no.  3,  Aug  1983. 

9.  The  Exciting  Mechanism  of  the  Parallel-Plate  EMP  Simulator,  H.M.  Shen, 
R.W.P.  King,  T.T.  Wu,  IEEE,  EHC-29,  no.  1,  Feb  1987. 

10.  Private  communication  FEL-TNO,  The  Haque,  The  Netherlands. 

11.  Sapiens  Memo  11,  Electromagnetic  characteristics  of  the  conical 
transmission-line  EMP  simulator  SAPIENS,  S.  Garmland,  FOA  Report  C 
30415-A3,  Feb  1986. 

12.  SSN  21,  Impedances  and  field  distributions  for  parallel  plate 
transmission  line  simulators.  C.2.  Baum,  AFWL,  6  June  1966. 

13.  SSN  52,  A  parameter  study  of  two  parallel  plate  transmission  line 
simulators  of  EMP  SSN  21,  T.L.  Brown,  K.D.  Granzow,  Dikewood 
Corporation,  19  April  1968. 

14.  Engineering  Techniques  for  Electromagnetic  Pulse  Hardness  Testing. 

A. L.  Whitson,  SRI,  Menlo  Park,  CA,  Sept  1974,  AD-786  722 

15.  Time-Domain  Analysis  of  a  Large  EMP  Simulator,  C.  Bardet,  0.  Dafif, 

B.  Jecko,  IEEE-EMC-29,  no.  1,  Feb  1987. 

16.  AFWL  Technical  Report  No.  AFWL-TR-69-80 


-  37  - 


17.  De  effecten  van  EMP  op  electronische  systemen,  P.A.A.  Sevat,  NERG, 
deel  48,  nr. 2,  1983. 

18.  Selleck,  B.,  AVRE,  Aldermaston,  UK,  Private  Communication,  1970. 

19.  Pulspak  8080,  A  transportable  all  weather  500  kV  EMP  generator,  Pulsar 
PATP  80-5-1,  Oct  1980. 

20.  Pulspaks  for  EMP,  Radiating  Pulse  Generator  Systems,  Phys.  Int.  Co, 

PTP  82-131,  Sept  1982. 

21.  I.A.D.  Lewis  and  F.H.  Veils,  Millimicrosecond  Pulse  Techniques, 
Pergamon  Press,  1959. 

22.  Radiation  from  idealized  shock  excitation  currents  in  a  straight 
conductor  rising  from  a  perfect  earth  at  an  arbitrary  angle, 

E.A.  Lewis,  Electromagnetic  Wave  Propagation,  Academic  Press,  1960. 


-  38  - 


APPENDIX  A 

ELECTROMAGNETIC  RADIATION 


A. 1  Introduction 


The  electromagnetic  radiation  from  parallel-plate  simulators  can 
roughly  be  approximated  by  the  following  method  (18).  When  a  step  function 
voltage  V  is  applied  at  the  source  end,  a  voltage  wave  of  amplitude  V  and 
current  V/Z^  starts  up  the  structure  vith  a  velocity  v,  close  to  the  velocity 
c  of  light.  The  propagation  of  this  transient  current  along  the  vires  sets  up 
radiation  fields. 


At  time  i/v  the  wave  reaches  the  horizontal  section,  at  2*i/v  the 
output  taper  and  at  3* i/v  the  termination  (see  Pig  Al).  The  radiation  of 
transients  by  straight  conductors  inclined  at  an  arbitrary  angle  a  from  the 
zenith  above  a  perfectly  conducting  ground  plane  has  been  described  by  Lewis 
(21,22).  Therefore  consider  first  the  input  taper  only. 


In  fig.  A2  the  conductor  is  in  the  plane  YOZ  and  the  observation  point 
P  has  the  spherical  co-ordinates  R,  0  and  ♦.  Levis  derived  the  following 
expressions  for  the  E  field  at  a  distance  R  from  the  source: 


E  _  _ 2kP*cos« 

E0  ”  1  -  0*sin+*sina 

„  -k.3>-sin(2g),cos8 

*  1  -  81 ■ cos J0* cos1 a 

Er-0 

where 

.  V 

AnCo’fe'Z^c 


(Al) 

(A2) 


(A3) 


0  is  the  ratio  of  v/c 
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In  table  A1  the  gain  factor  G  is  given  for  azimuth  angles  betveen  +90° 
(direction  of  current  flow),  through  0#  (broadside  to  simulator  axis)  to  -90°. 
Also  the  pulse  width  At  is  tabulated. 

The  product  G*At  in  (A7)  has  a  constant  value  i/v  -  10“^  s,  for  the 
worst  case  v  -  c.  In  practice  the  velocity  v  along  the  taper  is  less  than  c, 
so  At  in  (A7)  is  a  little  larger.  In  cases  where  At  is  smaller  than  the  rise 
time  tr  of  the  simulator  (pulser  and  taper)  only  a  part  of  the  amplitude  is 
built  up.  In  our  case  t  ■  5  ns.  A  correction  is  made  with  a  factor  At/5.  In 
figs.  A3  and  A4  the  field  patterns  are  plotted  for  the  6  m  and  10  m  simulator 
respectively,  and  in  figs.  A5  and  A6  the  corresponding  energy  density.  The 
dashed  lines  are  corrected  for  the  rise  time. 
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The  width  of  the  pulse  observed  at  a  distance  R  with  an  azimuth  angle 
♦  -  90®  is  very  small  because  the  velocity  of  the  current  in  the  taper  is 
assumed  to  be  the  same  as  the  velocity  of  the  propagation  of  the  radiated  wave 
and  because  the  difference  in  path  length  Ai  in  fig.  23  is  very  small.  In 
fig.  A7  the  timing  of  the  pulses  at  several  azimuth  angles  is  given. 

A. 3  Examination  of  the  Results 


The  radiation  levels  of  the  6  m  and  10  m  lines  were  calculated  by 
replacing  the  tapers  with  a  single  wire;  by  using  a  step  function  (t  =  0);  by 
assuming  a  perfectly  conducting  ground  and  no  propagation  loss  (Eg  *  1/R). 

In  practice  measurements  on  scale  models  and  on  real  EMP  simulators 
have  shown  that  the  levels  are  smaller.  At  azimuth  angles  +  -  -90°  the 
calculated  and  measured  peak  values  of  the  field  strength  are  about  the  same, 
at  =  0°  the  measured  values  may  be  a  factor  2-3  smaller,  and  at  ♦  90®  even 
more  than  a  factor  of  10.  That  is  because,  particularly  the  high  frequencies, 
are  more  than  linearly  attenuated  with  increasing  distance  R,  in  this  manner 
increasing  the  rise  time  and  reducing  the  radiation. 


A. 4  Radiation  Levels 


In  UK  and  the  Netherlands  the  following  limits  have  been  used  for  EHP 
Simulation  (18): 

-  Pace  makers,  400  V/m  repetitive,  1  pps 

-  EED's,  flash-over  1  kV/m, 

energy  12.5  uJ/ra*. 

The  minimum  distances  to  apply  for  safety  reason  can  be  assessed  by 
interpolating  or  extrapolating  the  figs.  A3  to  A6. 


-  42  - 


Figure  1:  Several  types  of  bounded-wave  simulators. 
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Wave  guide 


Transmission  line 


Reflections  by  local  discontinuities  of  Z\ 


Figure  2:  Tint  space  concept  for  single  pulse  excitation 
(no  SWR  but  clear  times). 
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Figure  3:  Model  for  exo-atmospheric  EMP. 
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Figure  4:  Spectrum  of  model  In  Figure  3. 


Figure  6:  Cumulative  energy  density  of  model  Figure  3. 


FIGURE  7:  FIELD  AND  POTENTIAL  DISTRIBUTION  FOR  SEMI -INFINITE,  PARALLEL, TWO-PLATE 
TRANSMISSION  LINE  (UPPER  HALF! 

Z  -  (W  +  1  +  ew)/  n  Z  =  x  +  jy  and  W  =  u  +  jv 


-  47  - 


^  2o 

) 

< 

| 

A.  INFINITE  OUTER  PLATES,  TI’REE  CONDUCTOR,  PARALLEL  PLATE 
TRANSMISSION  LINE 
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CONDUCTOR  PLANE  (iNfiNITE) 
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8.  TWO  CONDUCTOR,  PARALLEL  PLATE  TRANSMISSION  LINE 


FIGURE  8:  SYMMETRICAL,  PARALLEL  PLATE  TRANSMISSION  LINES 


FIGURE  9:  FIELD  PND  POTENTIAL  DISTRIBUTION 

PGR  PARALLEL .  TAO-PLATE  TRANSMISSION  LINE.  S7.97  OHMS 

3/P  =  .20 


FIGURE  10:  r  PELD  PND  POTENTIAL  DISTRIBUTION 

:'CR  PARALLEL.  TWO- PLATE  TRANSMISSION  LINE .  130-H9  OHMS 

B/A  =  .60 
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FIGURE  11:  r  Ic.L D  AND  POTENTIAL  DISTRIBUTION 

FOR  PARALLEL.  TWO-PLATE  TRANSMISSION  LINE.  1MM.03  OHMS 

B/A  =  .70 
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figure  13:  FIELD  AND  POTENTIAL  DISIRI3UTIGN 

FOR  PARALLEL.  TUG-PLATE  TRANSMISSION  LU'E.  167.71  OHMS 

B/A  =  .90 
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FIGURE  14:  FIELD  AND  POTENTIAL  DISTRIBUTION 

FOR  PARALLEL .  TWO-PLATE  TRANSMISSION  LINE.  178.18  OHMS 

B/A  =  1.00 
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FIGURE  15:  FIELD  AND  POTENTIAL  DISTRIBUTION 

FOR  PARALLEL.  TWO-PLATE  TRANSMISSION  LINE.  187.90  OHMS 

8/A  =  1 *10 
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FIGURE  16:  FIELD  AND  POTENTIAL  DISTRIBUTION 

FOR  PARALLEL.  TWO-PLATE  TRANSMISSION  LINE.  196.96 

B/A  =  1 .20 


figure  17:  FIELD  AND  POTENTIAL  DISTRIBUTION 

FOR  PARALLEL.  TWO-PLATE  TRANSMISSION  LINE.  205-H4  OHMS 

B/A  =  1 -30 
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FIGURE  18:  FIELD  AND  POTENTIAL  DISTRIBUTION 

FOR  PARALLEL.  TWO-PLATE  TRANSMISSION  LINE.  213. Ml  OHMS 

3/A  =  1  -MO 
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FIGURE  19:  FIELD  PND  POTENTIRL  DISTRIBUTION 

C0R  PARALLEL.  TWO-PLATE  TRANSMISSION  LINE.  253.02  OHMS 

B/A  =  2.00 
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FIGURE  20:  FIELD  AND  POTENTIAL  DISTRIBUTION 

FOR  PARALLEL.  TWO-PLATE  TRANSMISSION  LINE.  360. 08. OHMS 
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Figure  22:  Field  variations  across  a  two  plate  transmission  line. 


reflected  waves 


Figure  23:  Wave  front  field  structure  in  a  transmission  line  simulator. 


Figure  24:  Transmission  line  Impedance  as  a  function  of  the  lower 
plate  width. 


FIGURE  25: 


WIRE  GRIO  (ONE  CELL)  TERMINATING  UNIFORM 
FIELD  ON  ONE  SIDE  OF  GRID 


Figure 


d  enhancement  at  edge  of  semi-infinite  plate. 


Elevation  angle  0*  11.5° 
Apex  angle  4  •  1 8.6  0 


Figure  30:  Azimuth  and  elevation  angle. 


Correction  for  constant  impedance  of  the  launching  taper. 


167cm 
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INTERFACE 


Figure  32:  Dimensions  of  the  terminating  network. 


5100cm 
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Figure  35:  Example  of  a  wood  frame  work  (here  17  x  11  resistors). 


Figure  36:  Some  details  of  a  wooden  frame  work. 
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plates 


Figure  41:  Terminating  network  in  one  of  the  tapers 


Figure  42:  The  place  and  dimensions  of  the  terminating  networks  of 
the  10  m-line. 
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Figure  Als  The  parallel -plate  structure. 
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Figure  A2:  The 
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Figure  A3: 
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b/a  (dif) 

b/a  ^  (dif) 

b/a  f*  (dif) 

.010  0.009768 
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.070  0062229 

.000955 

.000870 
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.011  0.010723 
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.000867 

.000804 

.012  0.011675 
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.072  0.063838 

.000948 

.000865 

.000801 
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.000800 

.014  0.013569 
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.000941 

.000860 

.000798 

.015  0.014510 

.045  0.041430 

.075  0.06627 7 

.000939 

.000858 

.000796 

.016  0.915449 

.046  0.042288 

.076  0.067033 

.000935 

.000856 

.000794 

.017  0. 016384 

.047  0.043144 

.077  0.067827 

.000932 

.000853 

.000792 

.018  0.017316 

.048  0.043997 

.078  0.068G19 

.000929 

.000851 

.000791 

.019  0.018245 

.043  0.044848 

.079  0.069410 

.000926 

.000849 

.000789 

.020  0.019171 

.050  0.045697 

.080  0.070199 

.000923 

.000847 

.000786 

.021  0.020094 

.051  0.046544 

.081  0.070985 

.000920 

.000844 

.000786 

.022  0.121014 

.052  0.047388 

|  .082  0.071771 

.000918 

.000842 

.000733 

.023  0.021932 

.053  0.048230 

.083  0.072554 

.00097.4 

.000840 

.000781 

.024  0.022846 

.054  0.049070 

.084  0. 073^35 

.000911 

.000838 

.000780 

.025  0.023/57 

!  .055  0,049908 

i  .085  0.0/4115 

.000903 

.  .000336 

.000778 

.026  0.024666 

.056  0.050744 

.086  0.074893 

.000905 

.000834 

.000776 

.027  0.025571 

.057  0.051578 

.087  0.075669 

.000903 

-.000831 

.000775 

.028  0.026474 

.058  0.052409 

.088  0.076444 

.000901 

.000830 

.000773 

.029  0.027375 

.059  0.053239 

.089  0.077217 

.000897 

.000827 

.000771 

.030  0.028272 

.060  0.054066 

.090  0.0/7988 

.000395 

.000825 

.000769 

.031  0.029167 

.061  0.054891 

.091  0.078757 

.000893 

.000824 

. 000768 

.032  0.030060 

.062  0.055715 

.092  0.079525 

.000889 

.000821 

.000766 

.033  0.030949 

.063  0.056536 

.093  0.080291 

.000887 

.000819 

.000764 

.034  0.031836 

.064  0.057355 

.094  0.081055 

.000885 

.000817 

.000763 

.035  0.032721 

.065  0.058172 

.095  0.081818 

.000882 

.000313 

.000761 

.036  0.033603 

.066  0.058987 

.096  0.082579 

.000879 

.000814 

.000759 

.037  0.034482 

.067  0.059801 

.097  0.083338 

.000877 

.000811 

.000758 

.038  0.035359 

.068  0,060612 

.098  0.084C96 

.000875 

.000809 

.000756 

.039  0.036234 

.069  0.061421 

.099  0.084852 

.000872 

.000808 

.000758 

Table  1:  Geometric  factor  f  versus  b/a  =  0.01  -  0.1. 


b/» 

tl 

(dif  > 

0.10 

0.08561 

.00745 

0.11 

0.09306 

.00731 

0.12 

0.10037 

.00716 

0.  15 

0.10753 

.00702 

0.14 

0.11455 

.00689 

0.15 

0.12144 

.00676 

0.16 

0.12820 

.00663 

0.17 

0.13483 

.00653 

0.18 

0.14136 

.00641 

0.19 

0.14777 

.00630 

0.20 

0.15407 

.00619 

0.21 

0.16026 

.00610 

0.22 

0.16636 

.00600 

0.23 

0.17236 

.00591 

0.24 

0.17827 

.00581 

0.25 

0.10408 

.00573 

0.26 

0.18981 

.00504 

0.27 

0.19545 

.00556 

0.28 

0.20101 

.00548 

0.  29 

0.20649 

.00540 

0.30 

0.21189 

.00532 

0.31 

0.21721 

.00526 

0.32 

0.22247 

.00518 

0.33 

0.22765 

.00511 

0.34 

0.23276 

.00505 

0.35 

0.23781 

.00497 

0.36 

0.24278 

.00492 

0.37 

0.24770 

.00485 

0.38 

0.25255 

.00479 

0.39 

0.25734 

.00473 

b/a 

f 

Cdif ) 

0.40 

0. 26207 

.00468 

0.41 

0.2667  5 

.00462 

0.42 

0.27137 

.00456 

0.43 

0.27593 

.00451 

0.44 

0.28044 

.00445 

0.45 

0. 284S9 

.00441 

0.46 

0.28930 

.00435 

0.47 

0.  29365 

.00431 

0.48 

0.29796 

.00425 

0.49 

0.30221 

.00421 

0.  SO 

0.30642 

.00417 

0.51 

0.31059 

.00411 

0.52 

0.31470 

.00408 

0.53 

0.31878 

.00403 

0.54 

0.32281 

.00399 

0.55 

0.32680 

.00354 

0.56 

0.33074 

.00391 

0,57 

0.33465 

.00386 

0.58 

0.33851 

.00383 

0.59 

0.34234 

.00379 

0 .60 

0.34613 

.00375 

0.61 

0.34988 

.00371 

0.62 

0.35359 

.00368 

0.63 

9.35727 

.00364 

0.64 

0.36091 

.00361 

0.65 

0.36452 

.00357 

0.66 

0.36309 

.00354 

0.67 

0.37163 

.00350 

0.61 

0.37513 

.00347 

0.69 

0.37860 

.00344 

b/a 

Cdif) 

0.70 

0.38204 

.00341 

0.71 

0.38545 

.00338 

0.72 

0.38883 

.00335 

0.73 

0.39218 

.00331 

0.74 

0.39549 

.00329 

0.75 

0.39878 

.00326 

0.  76 

0.40204 

.00323 

0.77 

0.40527 

.00320 

C.  78 

0.40847 

.00318 

0.79 

0.41165 

.00314 

0.80 

0.41479 

.00313 

0.81 

0.41792 

.00309 

0.82 

0.42101 

.00307 

0.85 

0.42408 

.00304 

0.84 

0.42712 

.00302 

0.85 

C. 43014 

.00299 

0.83 

0.43313 

.00297 

0.87 

0.45610 

.00295 

0.88 

0.43905 

.00292 

0.89 

0.44197 

.00290 

0.90 

0.44487 

.00287 

0.91 

0.44774 

.00285 

0.92 

0.45059 

.00283 

0.93 

0.45342 

.00281 

0.94 

0.45623 

.00279 

0.95 

0.45902 

.00276 

0.96 

0.46178 

.00275 

0.97 

0.46453 

.00272 

0.98 

0.46725 

.00271 

0.99 

0.40996 

.00268 

Table  2:  Geometric  factor  f  versus  b/a  *  0.1  -  1.0. 
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b/a 

{_Z  (dif) 

1.00 

0.47264 

.01312 

1.05 

0.48576 

.01266 

1.10 

0.49842 

.01222 

1.15 

0. 51064 

.01181 

1.20 

0. 52245 

.01143 

1.25 

0. 53388 

.01107 

1.30 

0.54495 

.01073 

1.35 

0.55568 

.01041 

1.40 

0.56609 

.01010 

1.45 

0. 57619 

. C0382 

1. 50 

0.58601 

.00954 

1.55 

0.59555 

.00928 

1.60 

0.60483 

.00904 

1.65 

0.61387 

.00880 

1.70 

0.62267 

.00859 

1.75 

0.63126 

.00836 

1.80 

0.63962 

.C0817 

1.85 

0.64779 

.00797 

1.90 

0.65576 

.00779 

1.515 

0.66355 

.00761 

2.00 

0.67116 

.00745 

2.05 

0.67861 

.00728 

2.10 

0.68589 

.00712 

2.15 

0.69301 

.00698 

2.20 

0.59999 

.00683 

2.  25 

0.70682 

.00669 

2.30 

0. 71351 

.00656 

2.35 

0.72007 

.00643 

2.40 

0.72650 

.00632 

2.45 

0.73282 

.ooon 

Table  3 


b/a 

iz  (dif) 

2.50 

0.73901 

.00608 

2.55 

0.74509 

. 00596 

2.60 

0.75105 

. 00586 

2.65 

0.  75691 

.00576 

2.70 

0.76267 

.00566 

2.75 

0.76833 

.00556 

2.  80 

0.77389 

. 00547 

2.85 

0.77936 

. 00538 

2.  90 

0.78474 

. 00530 

2.9S 

0.79004 

. 00S21 

3.00 

0.79525 

.00512 

3.05 

0.80037 

.00505 

3.10 

0.80542 

. 00497 

3. ’5 

0.81039 

.00490 

3.  20 

0.81529 

.00482 

3.25 

0.82011 

. 00475 

3.  30 

0.82486 

.00469 

3.35 

0.82955 

.00461 

3.  40 

0.83416 

.00456 

3.45 

0.83872 

.00449 

3.50 

0.84321 

.00443 

3.55 

0.84764 

.00436 

3.60 

0.85200 

.00431 

3.65 

0.85631 

.00426 

3.70 

0.86057 

.00420 

3.75 

0.86477 

.00414 

3.80 

0.86891 

.00410 

3.85 

0.87301 

.00404 

3.90 

0.87705 

.00399 

3.95 

0.88104 

.00394 

Geometric  factor  f 

9 


b/a 

f_Z  (dif) 

4.00 

0.88498 

.00390 

4.05 

0.88888 

.00385 

4.10 

0.89273 

.00380 

4  .15 

0.89653 

.00376 

4.20 

0.90029 

.00372 

4.25 

0.90401 

.00367 

4.30 

0.90768 

.00363 

4.35 

0.91131 

.00359 

4.40 

0.91490 

.00356 

4.45 

0.91846 

.00351 

4.  50 

0.92197 

.00347 

4.55 

0,92544 

.00344 

4.60 

0.92888 

.09341 

4.65 

0.93229 

.00336 

4.70 

0.93565 

.  0  "  3  3 

9.75 

0.93893 

.00330 

4.80 

0.94228 

.00327 

4.85 

0.94555 

.00323 

4.90 

0.94878 

.00320 

4.95 

0.95198 

.00316 

5.00 

0.95514 

.00314 

5.05 

0. 95828 

.00311 

5.  1C 

0.96139 

.00307 

5.15 

0.96446 

.00305 

5.  26 

0. 96751 

.00302 

5.25 

0.97053 

.00299 

5.30 

0.97352 

.00296 

5.35 

0.97648 

. 09294 

5 .  4  C 

0.9794’ 

. CC20 1 

6.4  4 

0.98233 

.00287 

sus  b/a  *  1.0  -  5.5. 

/ 


-  90  - 

b/a  *g  (dif) 

b/a  fR  (dif) 

b/a  fg  (dif) 

5.50  0.9852 

7.00  1.0615 

8.50  1.1230 

.  0029 

,  .0022 

. . .  .0019 

5.55  0.9881 

7.05  1.0637 

8.55  1.1249 

.0028 

.  0023 

.0018 

5.60  0.9909 

7.10  1.0660 

8.60  1.1267  ▼ 

.  0028 

.0022 

.0019 

5.65  0.9937 

7.15  1.0682 

8.65  1.1286 

.  0028 

.  0022 

.0018 

5.70  0.9905 

7.20  1.0704 

8.70  1.1304 

.  0028 

.  0022 

.0018 

5.75  0.9993 

7.25  1.0726 

8.75  1.1322 

.0027 

.0022 

.0018 

5.80  1.0020 

7.30  1.0748 

8.80  1.1340 

.0027 

.  0021 

.0018 

5.85  1.0047 

7.35  1.0769 

8.85  1.1358 

.0027 

.0022 

.0018 

5.90  1.0074 

7.40  7.0791 

8.90  1.1376 

.0027 

.  0021 

.0018 

5.95  1.0101 

7.45  1.0812 

8.95  1.1394 

.0026 

.0021 

.0018 

6.00  1.0127 

7.50  1.0833 

9.00  1.1412 

.  0026 

.0021 

.0017 

6.05  1.0153 

7.55  1.0854 

9.05  1.1429 

.0026 

.0021 

.0018 

6.10  1.0179 

7.60  1.0875 

9.10  3.1447 

.  0026 

.0021 

.0017 

6.15  1.0205 

7.65  1.0896 

9.15  1.1464 

1  .0026 

.  0021 

.0017 

6.20  1.0231 

7.70  1.0917 

9.20  1.148i  g 

.  0C25 

.0020 

.0018  & 

6.25  1.0256 

7,75  1.0937 

9.25  1.1499  S 

.0025 

.  0021 

. OC  i  7  '  K 

6.30  1.0281 

7.80  1.0958 

9.30  1.1516 

.  0025 

.0020 

.0017 

6.35  1.0306 

7.85  1.09  78 

9.35  1.1533 

.  0025 

.0020 

.0017 

6.40  1.0331 

7.90  1.0998 

9.40  1.1550 

.  0025 

.0020 

.0016 

6.45  1.035G 

7.95  1 . 1C  18 

9.45  1.1566 

.0024 

.  0020 

.0017 

6.50  1.0380 

8.00  1.1038 

9.50  1.1583 

.0024 

.  0020 

.0017 

6.55  1.0404 

8,05  1.1058 

9.55  1.1600 

.  0025 

.0019 

.0010 

6.60  1.0429 

8.10  1.1077 

9.60  1.1616 

.  0023 

.  0020 

.0017 

6.G5  1.0452 

8.15  1.1097 

".65  1.1633 

.0024 

.0019 

.0016 

6.70  1.0476 

8.20  1.1116 

9.70  1.1649 

.  0024 

.0020 

.0017 

6.75  1.0500 

8.25  1.1136 

9.75  1.1666 

.0023 

.0019 

.0016 

6.80  1.0523 

8.30  1.1155 

9.80  1.1682 

.0023 

.0019 

.0016 

6.85  1.0546 

8.35  1.1174 

9.85  1.1698 

.  0023 

.0019 

.0016 

6.90  1.0569 

8.40  1.1193 

9.90  1.1714 

.  0023 

.0019 

.0016 

6.93  1.0592 

8.45  1.1212 

9.95  1.1730 

.0023 

.0018 

.0016 

Table  4:  Geometric  factor  f  versus 

9 

b/a  »  5.5  -  10 

\ 

\ 
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(dif) 


(dif) 
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TABLE  6 


h 

(m) 

• 

r/h 

f 

Ef 

(kV/cn) 

E 

e 

(kV/cn) 

6 

6.7 

I 

o 

H 

22 

0.75 

16.5 

5 

8 

»-* 

o 

1 

4S 

20 

0.9 

18 

4 

1 

1 

O 

H 

18 

1.13 

20.3 

3 

1.3 

»-» 

o 

1 

U> 

15.5 

1.5 

23.3 

2 

2 

•  10-3 

12.5 

2.25 

28.1 

1 

4 

•  10"3 

9 

4.5 

40.5 

Table  6:  The  edge  enhanced  field  Ee  of  the  taper 
(wire  radius  r  *  4  mm). 


Table  7.  Enhanced  field  at  the  edge  Eg  and  between  the  vires 
(edge  radius  is  5  cm,  wire  radius  is  4  mm). 
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Table  9.  Enhanced  field  E  at  the  edge  and  E  between  the  vires 

«?  w 

of  the  lOm-line. 


TABLE  10 
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Table  10.  Correction  of  height  of  the  taper  for  constant 
of  the  lOa-line. 
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Table  11 


Specifications 

6m-llne 

lOm-line 

dimensions  l  x  v  x  1 

i  90  x  10  x  6  m 

90  x  20  x  10  m 

working  volume 

(15-30)  x  10  x  3  m3 

(15-30)  x  20  x  5  m 

ground  plate 

vires,  v  »  20  m 

wire  mesh,  20  m 

line  impedance 

101.5  a  +  14  a 

90.7  a  +  24  2 

rise  time 

(  elevation  2  ns 

5.4  ns 

1 

{  azimuth  1.4  ns 

5.4  ns 

field  magnitude 

75  kV/m 

25-50  kV/m 

field  distortion  'j 

|  502 

502 

field  variation  J 

field  uniformity  'J 

|  fig.  16 

fig.  14 

field  orientation  J 

Table  11.  Summary  of  the  specifications  of  the  6m  and  10m  line. 
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TABLE  A1 


angle  4  (°) 

6m- 

G 

■  line 

At(ns) 

10m- 

G 

•line 

At(ns) 

-90 

0.5 

198 

0.51 

195 

-75 

0.522 

192 

-60 

0.54 

185 

0.55 

182 

-45 

0.6 

167 

-30 

0.67 

149 

0.68 

147 

-15 

0.8 

125 

0 

1 

100 

1 

100 

15 

1.34 

74.6 

1.33 

75 

30 

1.96 

51 

1.9 

53 

50 

4.02 

24.9 

3.66 

27.3 

60 

6.63 

15.1 

5.6 

17.8 

65 

8.99 

11.1 

70 

12.7 

7.9 

9.23 

10.8 

75 

18.9 

5.3 

12 

8.34 

80 

19.6 

3.4 

15.3 

6.6 

85 

18.3 

5.47 

90 

19.6 

1.9 

19.6 

5.1 

Table  Al.  The  gain  factor  G  and  pulse  width  At  for  the 
6m  and  10a  line  with  rise  time  tf  »  2  ns  and 
5.4  ns  respectively 
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